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CMB power spectra and the SNIa luminosity-redshift curves have difficulty distinguishing between 
models with the same average value of the dark energy equation of state. We propose using the 
cross-correlation of the CMB temperature anisotropy with large scale structure to help break this 
degeneracy. 
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The evidence for Dark Energy (DE), a dark compo- 
nent causing the accelerated expansion of the universe, 
comes from several complementary sources. The analysis 
of the cosmic microwave background (CMB) shows that 
the total energy density of the universe is very close to 
the critical density, i.e. the universe it flat. At the same 
time, large scale structure (LSS) measurements suggest 
that no more than a third of the critical energy density 
can be in the form of clustering matter. While CMB 
and LSS measurements, if considered separately, can be 
explained without invoking DE, they require DE for con- 
sistency with each other. 0,01 • A more direct evidence 
for DE comes from measurements of the luminositydis- 
tance vs redshift for supernovae type la (SNIa) 0, Q 
revealing an accelerated expansion of the universe. For a 
flat Freedman-Robertson- Walker (FRW) universe accel- 
eration implies domination by a substance with negative 
pressure, such as vacuum energy. Another direct evi- 
dence which has recently become available is the detec- 
tion of the Integrated Sachs- Wolfe (ISW) effect using the 
CMB/LSS cross-correlation H H H @, 13 E3 ■ This 
evidence is of complementary nature to the SNIa data. 
Rather than probing the accelerated expansion it reflects 
the slow down in the growth of density perturbations ex- 
pected in a non-matter dominated universe. 

The cosmological constant A is the simplest DE model 
giving a satisfactory fit to the existing data. Observa- 
tions are also consistent with an evolving dark energy, 
such as the scalar field Quintessence 0, [lJ 0] . Estab- 
lishing whether the dark energy is constant or evolving 
is one of the main challenges for modern cosmology. In 
the FRW universe the large scale evolution of DE is de- 
termined by its equation of state defined as the ratio of 
pressure to energy density: w = p/p. For scalar field 
Quintessence w also determines the clustering properties 
of DE. Depending on the model w can be constant or 
change with time. Models with w ^ — 1 correspond to 
evolving DE, while w = — 1 corresponds to A. 

Much work has been done on trying to constrain the 
change in w by fitting various forms of w(z) to the 
SNIa luminosity distance-redshift data, often in com- 
bination with constraints from the CMB power spec- 
trum. It is known, however, that these observables de- 
pend on w(z) via one or more integrals over the red- 
shift [HIIlIilEEIi^. As 

we illustrate below, the 
CMB power spectra and the luminosity distance curves 



for models with varying w(z) can be difficult to distin- 
guish from predictions of the model with a constant w 
equal to the average of w(z) defined as 



< w >z 
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where a is the scale factor, a = 1 today, a* is the scale 
factor at the last scattering and Clrt(a) is the ratio of the 
DE energy density to the critical density. 

The LSS measurements probe the net growth of cosmic 
structures. Provided other parameters are known, LSS 
can differentiate between a constant and a time-varying 
w, especially if one had the information about the growth 
factor at different redshifts |22l |. 

Measurements of the ISW effect can provide another 
probe of the evolution of w(z). In a way, they are more 
sensitive to the change in w(z) because ISW depends on 
the growth rate of cosmic structures (the time-derivative 
of the growth factor) in addition to the total growth (see 
[23| for a useful discussion). In this letter we build on 
the methods and results of [24| and illustrate the poten- 
tial of the CMB/LSS cross-correlation for constraining 
the change in w(z). For this purpose we will consider 
four Quintessence models with the same average < w > 
shown in Fig. ^ These four models are deliberately cho- 
sen to illustrate a point and do not necessarily represent 
particularly well-motivated quintessence potentials. The 
ansatz for w(z) used to construct our models is: 
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It is similar in spirit to that of |25j and describes a tran- 
sition from w + = w{z = +00) to = w{z = —00) with 
parameters Sz and zq describing the width and the cen- 
tral redshift of the transition. The values of the ansatz 
parameters used for the four models in Fig.0are given in 
Table [I] The shapes of the corresponding Quintessence 
potentials can be reconstructed from w{z), if it was de- 
sired. These four models have one thing in common - 
they all have the same < w >— —0.8, as defined in 
eq. |0. This value, < w >— —0.8, corresponds to the 
upper boundary of the 95% confidence region obtained 
from the combined analysis of the WMAP, SDSS and 
SNIa data (Fig. 13 of 0). 
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dot-dash 
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TABLE I: Models considered in the paper. 



FIG. 1: w vs z for the four varying w(z) models considered 
in this paper. They all have the same average value < w > 
(solid line) given by cq.0 



In Fig. [3 we show the CMB spectra for the models in 
Table [|] The cosmological parameters are the same for 
all models. Through the paper we assume a flat universe 
with the Hubble parameter h — 0.66, baryon density 
n b h 2 — 0.024, total matter density VL M h 2 = 0.14, spec- 
tral index n s — 0.99, optical depth r = 0.166 and the 
amplitude of scalar fluctuations A s — 0.86 (as defined in 
|26|.) As evident from Fig. [21 all four varying w models fit 
the WMAP data very well and are indistinguishable from 
each other and from the model with a constant w = —0.8. 
This agrees with other studies (e. g. [2 ^) an d supports 
the point made in (among other papers) [lg that CMB 
spectra are mainly sensitive to the averaged value of w, 
as defined in eq. 1 - 



1 The choice of the parameterization of w(z) has a major influence 
on the strength of the constraints. E. g., mod els with w(z) chosen 
to be linear in a or z disfavor rapid variation |28| . However, given 
sufficient freedom in a mode L curren t data is perfectly consistent 
with a rapidly varying w(z) l29tl3d . 




FIG. 2: CMB temperature angular power spectrum (TT) and 
the temperature-polarization cross-correlation (TE) for the 
five models in Fig. Q (using the same conventions for the dif- 
ferent line types) with the WMAP's first year data. 




FIG. 3: The percent difference in <1l(z) for the five models in 
Fig. [T]and the ACDM model (thick solid line). 



The degeneracy between the models persists to a large 
extent in the luminosity-redshift curves. The SNIa mea- 
surements probe the luminosity distance Hq(1l(z), where 



d L {z) = (l + z) 



z dz' 



(3) 



la H(z<) ' 

H (z) is the Hubble parameter and Hq is its current value. 
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FIG. 4: The linear matter power spectrum for the five models 
in Fig. H] The SDSS data points [53 are plotted to show 
current experimental error bars. 



In Fig. we plot the percent difference between each of 
the models with varying w and the model of constant 
W = —0.8. As can be seen from the figure, the differences 
from the w = const case do not exceed 5%, which is just 
about, if not slightly below, the current level of accuracy 
with which g?l(z) is extracted from SNIa. Therefore, a 
significantly varying w(z) can be consistent with SNIa if 
the model has a reasonable average < w >. Most opti- 
mistically, future SNIa data can improve the accuracy in 
determination of (Il(z) to the level of 1 — 2% ^7j- As 
the solid, dash-dot and long dash lines in Fig. |3] show, 
the future SNIa data may be able to distinguish between 
an increasing and a decreasing w(z), but there is still 
likely to be a constant w that gives a very similar for 
the same cosmological parameters. This point was pre- 
viously emphasized in 0,0. Fig. also shows a clear 
difference in predictions of the w = —0.8 model and the 
ACDM model. The ACDM model requires h = 0.72 for 
the same Ojf/i 2 in order to fit the CMB power spectra, 
hence the ACDM curve in Fig.|31corresponds to a smaller 

Fig. 0] shows the matter power spectra for the mod- 
els considered in the paper. Plotted are the linear CDM 
spectra at z = 0. A proper accounting for the non-linear 
effects as well as for the redshift space distortion intro- 
duces non-negligible modifications to the amplitudes and 
the shapes of the spectra. We did not attempt to do 
that here, assuming that changes in the differences be- 
tween the curves would be small. The SDSS data points 
are only plotted to show the uncertainty in the current 
determination of P(k). Not surprisingly, the models in 
which w decreases with time are easier to distinguish, be- 



cause the dark energy component becomes significant and 
starts to affect the growth of fluctuations early. On the 
other hand, in models with increasing w(z) the growth 
of fluctuations is not significantly affected until very re- 
cently and the change in w with time is harder to detect. 
Knowing the matter power spectrum at different redshifts 
could provide a useful handle on the time- variation of w 

M 

If one assumes a flat universe with adiabatic initial 
conditions then the combination of the CMB with SNIa 
data provides the information about the cosmological pa- 
rameters and < w >. The LSS data tightens the range of 
parameter values, fixes the galaxy bias and can give some 
information about the time variation of w. Measure- 
ments of the CMB/LSS cross-correlation can play a com- 
plementary role and provide further constraints on the 
time variation of w. What makes the cross-correlation 
interesting for this purpose is that it depends not only on 
the total growth of cosmic structure but on the growth 
rate as well. 

Cross-correlating the CMB with LSS with the purpose 
of detecting the ISW effect was first proposed in [2l( . Let 
us define 

A(n)=T(n)-f (4) 

and 

S(n) = , (5) 

where T(n) is the CMB temperature measured along the 
direction n, p(n) is the mass density along n, and T and 
p are the averaged CMB temperature and the matter 
density. The temperature anisotropy due to the ISW 
effect is an integral over the conformal time: 

-U-= J dr)^-*)^-^}, (6) 

where rj r is some initial time deep in the radiation era, 
i]o is the time today, $ and \& are the Newtonian gauge 
gravitational potentials, and the dot denotes differentia- 
tion with respect to rj. The quantity <5(n) contains con- 
tributions from astrophysical objects (e. g. galaxies) at 
different redshifts and can also be expressed as an inte- 
gral over the conformal time: 

f Vo dz 

S(n) = J dn- W g (z(v)) *(fa> - ri)n, v), (?) 

where W g {z) is a normalized galaxy selection function. 
We do not explicitly include a possible bias factor, assum- 
ing it is determined from the measurements of the LSS. 
We are interested in calculating the cross-correlation 
function 

X{6) = A(|nx - n 2 |) = (A(n 1 )5(n 2 )) , (8) 

where the angular brackets denote ensemble averaging 
and is the angle between directions fii and fi2. It can 
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be written as 124 



X{6) 



9 

25 



dm 



7] r J-q T 

dk . 9 N sm(kR) 



kR 



(9) 



where is the wave-number, A^(fc) is the primordial cur- 
vature power spectrum, F(k,r}\,r}2) describes the time- 
evolution of ($ — \j Q a nd 5, R = \Jr\ + r\ — 2rir 2 cos9 
and r = rjQ — rj (see [24| for more details) 2 . 

The galaxy window function W g {z) plays an important 
role as it essentially zooms in on the ISW effect over 
a particular redshift range. To help one decide on the 
selection of W g (z) it is useful to define 1(9, z) as 



1(9, z) - I d m I ^A^(fc) Fik . in . i/2 1 



_9_ '" /!! 
25.,, 



kR 



so that eq. J5J can be written as 
rO 



X(9) 



dz W g {z)I{z,0) 



(10) 
(11) 



Typically, the theoretical prediction for the cross- 
correlation approaches a flat plateau for 9 < 0.2° and 
the hight of the plateau is a good measure of the over- 
all signal. In Fig. we plot I(z,9 = 0.2°) vs z for the 
models considered in the paper. From this plot one can 
guess where to place the window function to see the max- 
imum difference in X(9) for the models. This, however, 
may not correspond to the choice that gives the high- 
est signal-to-noise ratio. Due to the arbitrary nature of 
the models considered in the paper we did not search 
for the optimal window function(s) that would maximize 
the difference between the models while minimizing the 
statistical uncertainty in X(9). Such optimization proce- 
dure, however, would be necessary if the cross-correlation 
was used to study properties of dark energy. 

For computational purposes it is advantageous to de- 
compose X(9) into a Legendre series, 



X(9) 



E- 

1=2 



1 



4tt 



-X e P e (cos 9) 



(12) 



In Fig. El we show results for the angular cross-correlation 
1(1 + 1)Xi/2tt for the two choices of Gaussian window 
functions shown in Fig. |SJ One is centered at z w = 0.3 
and the other at z w = 0.6, both have a standard de- 
viations of Az = 0.15. The two window functions are 



2 In addition, one has to subtract the monopole and dipole contri- 
butions to cq The details of how it was done can be found 
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FIG. 5: The integrand I(z,6 = 0.2°) for the five models in 
Fig. ^] Plotted on top (long dash-dot) are the two choices of 
window functions used in the paper, rescaled to fit the box. 



chosen so that the first one results in the larger differ- 
ence between the models, while the second one gives the 
larger signal-to-noise ratio. Theoretical uncertainties in 
individual Xg are too large to allow comparing the mod- 
els at each I. Instead one could calculate the total signal, 
Xt, for each model by summing (21 + l)Xt/4w over I: 



Xt — — : X(. 



i=i 



4tt 



(13) 



Note that, for Z max = oo, X T = X(9 = 0). The best 
achievable signal-to-noise in determination of Xt is ap- 
proximately given by 



sky / , 

(=2 



(2i + l) 



X] 



r'TTr'MM 



(14) 



where / s k y is the surveyed fraction of the sky, Cj T is the 
angular CMB temperature power spectrum and C^ M 
is the analogously defined matter angular spectrum in- 
side a given window W g (z). Table UTI contains Xt for 
each model for the two choices of the window functions. 
The table shows differences of up to ~ 35% in Xt for 
the window function centered at z w = 0.3 and up to 
~ 25% for the one centered at z w — 0.6. The signal- 
to-noise does not significantly depend on the particular 
model. Typically, S/N rs 2 for the window function with 
z w = 0.3 and S/N ph 3 for z w = 0.6, assuming complete 
sky coverage (/ s ky = 1) and that Cj T contains not just 
the ISW contribution but the entire CMB anisotropy on 
relevant scales. Neither choice of the window function 
would distinguish between the models considered in the 



5 



0.06 



t 0.04 

w 







z w =0.3 , Az = 0.15 . 


z w = 


0.6 , Az = 0.15 . 


//V 


// 


_ 

A 


-ill ! % 


/ 


i 

^\ 


/// 4 


■ 1 

II 


n 


\ 

1 1 1^1 


- #/ 

f 




«. 

\\ 

\- 


10 100 


10 


100 



1 1 

FIG. 6: The angular cross-correlation, 1(1 + 1)Xi/2tt for the 
five models in Fig. [I] and two choices of window functions. 
Theoretical uncertainties in individual Xt, if plotted, would 
be well outside the boundaries of the box. Instead, we pro- 
pose using the cumulative cross-correlation (summed over I) 
to differentiate between the models. 



Model/line type 


X T (z w = 0.3) 


X T (z w = 0.6) 


id — —0.8, thick solid 


0.167 


0.125 


dot 


0.187 


0.137 


short dash 


0.132 


0.108 


long dash 


0.175 


0.132 


dot-dash 


0.144 


0.106 



TABLE II: The cumulative cross-correlation Xt, in fiK, for 
< ma i = 300 for the models considered in the paper and the 
two choices of window functions. 



paper. There is hope, however, that the uncertainty in 
Xt can be reduced to the level that would allow model 
comparison. 

The potential accuracy that can be achieved in cross- 
correlation measurements was studied in The 
main theoretical uncertainty arises from limitations im- 
posed by cosmic variance and the fact that, in addition 
to the ISW contribution, the CMB signal has a sizable 
primary component. A possible way to increase the sig- 
nal to noise, proposed in |33|. is to consider the cross- 
correlation in multiple redshift shells. If the shells could 
be selected in a way that they were nearly uncorrelated 
then one could treat the cross-correlation in each shell 
as a separate measurement. Based on the results of [33] 
it appears that the best accuracy one could hope for is 
about 10%. If this was the case, the cross-correlation 
could be useful for adding to the constraints on rapidly 
varying w(z). The CMB/LSS correlation could also help 
with constraining the clustering properties of Dark En- 
ergy (see e.g. [34tl35| ). 

There is still room for a better understanding of the 
inherent limitations of the cross-correlation analysis. In 
particular, there is need for a procedure that would opti- 
mize the selection of redshift shells to maximize the signal 
to noise. This, as well as the constraints from the exist- 
ing data on the variation of w(z), is a subject of ongoing 
work |3(j. 
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